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1J.  AfSTRACT  (MsMunum  200  wofW 

This  Final  Report  describes  our  research  of  hypersonic  and  supersonic  boundary-layer 

flows. 

In  spite  of  its  extreme  importance  to  the  accurate  prediction  of  drag  and  heating 
requirements  in  high-speed  flow,  the  study  of  boundary- layer  transition  in  hypersonic  (NASP) 
and  supersonic  (fighter  and  high-speed  civil  transport)  flows  is  still  very  much  in  its  infancy. 
Transition  is  well  known,  however,  to  depend  strongly  on  such  effects  as  pressure  gradient,  wall 
curvature,  sweep,  roughness,  wall  mass  transfer,  freestream  and  wall  temperature,  nose  radius, 
nonequilibrium  chemistry,  and  freestream  disturbances.  (These  effects  have  been  discussed  in 
any  number  of  workshops  and  U.S.  Transition  Study  Group  meetings  under  the  direction  of  Eli 
Reshotko.) 

We  have  completed  detailed  studies  of  the  stability  of  the  laminar  basic  state  of  2-D  and 
axisymmetric  boundary  layers  with  non-equilibrium  chemistry  included  and  3-D  boundary-layer 
flows  of  an  ideal  gas.  (Relatively  simple  geometries  were  considered  due  to  the  anticipated 
difficulties  in  performing  basic-state  analyses.)  ~ 
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ABSTRACT 


This  Final  Report  describes  our  research  of  hypersonic  and  supersonic  boundary-layer 

flows. 

In  spite  of  its  extreme  importance  to  the  accurate  prediction  of  drag  and  heating 
requirements  in  high-speed  flow,  the  study  of  boundary-layer  transition  in  hypersonic  (NASP) 
and  supersonic  (fighter  and  high-speed  civil  transport)  flows  is  still  very  much  in  its  infancy. 
Transition  is  well  known,  however,  to  depend  strongly  on  such  effects  as  pressure  gradient,  wall 
curvature,  sweep,  roughness,  wall  mass  transfer,  freestream  and  wall  temperature,  nose  radius, 
nonequilibrium  chemistry,  and  freestream  disturbances.  (These  effects  have  been  discussed  in 
any  number  of  workshops  and  U.S.  Transition  Study  Group  meetings  under  the  direction  of  Eli 
Reshotko.) 

We  have  completed  detailed  studies  of  the  stability  of  the  laminar  basic  state  of  2-D  and 
axisymmetric  boundary  layers  with  non-equilibrium  chemistry  included  and  3-D  boundary-layer 
flows  of  an  ideal  gas.  (Relatively  simple  geometries  were  considered  due  to  the  anticipated 
difficulties  in  performing  basic-state  analyses.)  The  complete  region  between  the  wall  and  the 
shock  was  considered  in  the  non-equilibrium  chemistry  work. 

In  2-D  and  axisymmetric  flows,  the  inclusion  of  chemistry  caused  a  shift  of  the  second 
mode  of  Mack  to  lower  frequencies.  This  was  found  to  be  due  to  the  increase  in  size  of  the 
region  of  relative  supersonic  flow  because  of  the  lower  speeds  of  sound  in  the  relatively  cooler 
boundary  layers.  It  was  also  found  that  equilibrium  and  nonequilibrium  solutions  could  be  very 
different  depending  on  the  rates  of  the  reactions  relative  to  the  time  scales  of  convection  and 
diffusion.  In  particular,  in  equilibrium-air  calculations,  modes  which  travel  supersonically 
relative  to  the  inviscid  region  were  shown  to  exist.  These  modes  were  a  superposition  of 
incoming  and  outgoing  disturbances  which  exhibited  oscillatory  behavior  because  of  the  finite 
shock  stand-off  distance. 

For  the  examination  of  3-D  effects,  a  rotating  cone  was  used  as  a  model  of  a  swept  wing. 
An  increase  of  stagnation  temperature  was  found  to  be  only  slightly  stabilizing.  The  correlation 
of  transition  location  (N  =  9)  with  parameters  describing  the  crossflow  profile  was  investigated 
on  the  rotating  cone.  Transition  location  did  not  correlate  with  the  traditional  crossflow 
Reynolds  number.  A  new  parameter  that  appears  to  correlate  for  boundary-layer  flow  was 
found.  A  verification  with  the  NASA/Langley  Mach  3.5  Quiet  Tunnel  experiments  on  the 
yawed  cone  was  done. 
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1.  INTRODUCTION 

In  this  Final  Report,  Section  2  documents  our  technical  accomplishments.  Sections  3-9 
contain  more  detailed  information  of  the  results  of  our  research.  Section  10  describes  the 
resources  available  for  the  research,  and  Appendix  A  contains  a  short  resume  for  the  Principal 
Investigator. 

2.  TECHNICAL  ACCOMPLISHMENTS 

Students:  1  PostDoctoral  student,  2  PhD  students,  5  MS  students,  5  BS  students  involved  in 
related  work 

Journal  articles:  2  published  in  Annual  Review  of  Fluid  Mechanics,  2  published  in  Physics  of 

Fluids  A,  2  submitted  to  AIAA  Journal,  1  submitted  to  Journal  of  Computational  Physics 

Conference  papers:  2  invited,  6  other 

Conference  and  industry  presentations:  16  invited,  9  other 

Publications 

1.  "Stability  of  Three-Dimensional  Boundary  Layers,"  H.L.  Reed  and  W.S.  Saric,  Ann. 
Rev.  Fluid  Mech.,  Vol.  21,  p.  235,  Jan.  1989. 

2.  "Compressible  Boundary-Layer  Stability  Theory,"  H.L.  Reed  and  P.  Balakumar,  Phys. 
Fluids  A,  Vol.  2,  No.  8,  p.  1341,  Aug.  1990. 

3.  "Stability  of  Three-Dimensional  Supersonic  Boundary  Layers,”  P.  Balakumar  and  H.L. 
Reed,  Phys.  Fluids  A,  Vol.  3,  No.  4,  p.  617,  Apr.  1991. 

4.  "A  Catalogue  of  Linear  Stability  Theory  Results,"  H.L.  Reed,  accepted  Ann.  Rev.  Fluid 
Mech.,  Vol.  26,  1994. 

5.  "On  the  Linear  Stability  of  Supersonic  Cone  Boundary  Layers,"  G.K.  Stuckert  and  H.L. 
Reed,  accepted  AIAA  J,  1992. 

6.  "Linear  Disturbances  in  Hypersonic,  Chemically  Reacting  Shock  Layers,"  G.K. 
Stuckert  and  H.L.  Reed,  submitted  to  J  Comp.  Phys.,  1992. 

7.  "Transition  Correlation  in  Three-Dimensional  Supersonic  Boundary  Layers,"  H.L.  Reed 
and  T.S.  Haynes,  in  preparation  to  be  submitted  to  AIAA  J,  1992. 

8.  "Stability  and  Transition  of  Three-Dimensional  Boundary  Layers,"  W.S.  Saric  and  H.L. 
Reed,  Invited  Paper,  AGARD  Conference  No.  438,  Fluid  Dynamics  of  Three- 
Dimensional  Turbulent  Shear  Flows  and  Transition,  Cesme,  Turkey,  Oct.  1988. 

9.  "Supersonic/Hypersonic  Laminar A’urbulent  Transition,"  H.L.  Reed,  G.K.  Stuckert,  and 
P.  Balakumar,  Invited  Paper,  in  Developments  in  Mechanics,  Vol.  15,  Proceedings  of 
the  21st  Midwestern  Mechanics  Conference,  Aug.  13-16,  1989. 

10.  "Stability  of  High-Speed  Chemically  Reacting  and  Three-Dimensional  Boundary 
Layers,"  H.L.  Reed,  G.K.  Stuckert,  and  P.  Balakumar,  3rd  lUTAM  Symposium  on 
Laminar-Turbulent  Transition,  ed.  R.  Michel  and  D.  Amal,  Springer-Verlag,  New 
York,  1991. 

11.  "Stability  Limits  of  Supersonic  Three-Dimensional  Boundary  Layers,"  H.L.  Reed,  T. 
Haynes,  and  P.  Balakumar,  AIAA  Paper  90-1528. 

12.  "Stability  of  Hypersonic,  Chemically  Reacting  Viscous  Flows,”  G.K.  Stuckert  and  H.L. 
Reed,  AIAA  Paper  90-1529. 

13.  "Unstable  Branches  of  a  Hypersonic,  Chemically  Reacting  Boundary  Layer,"  G.K. 
Stuckert  and  H.L.  Reed,  Royal  Aeronautical  Society  Conference  on  Boundary-Layer 
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Transition  and  Control,  Cambridge  UK,  April  1991. 

14.  "Observations  in  Using  Linear  Stability  Theory  for  3-D  Supersonic  Boundary  Layers," 
H.L.  Reed  and  T.S.  Haynes,  Eli  Reshotko's  60th  Birthday  Party,  ICASE  and 
NASA/Langley  Research  Center,  July  28,  1991,  to  appear  1992. 

15.  "Stability  of  Hypersonic  Boundary-Layer  Flows  with  Chemistry,"  H.L.  Reed,  G.K. 
Stuckert,  and  T.S.  Haynes,  70th  AGARD  Fluid  Dynamics  Panel  Symposium  on 
Theoretical  and  Experimental  Methods  in  Hypersonic  Flows,  Torino,  Italy,  May  4-8, 
1992. 

Presentations 

1.  "Three-Dimensional  Boundary-Layer  Stability,"  H.L.  Reed,  Invited  Talk,  Naval  Post 
Graduate  School,  Monterey,  February  18,  1988. 

2.  "Computational  Simulation  of  Three-Dimensional  Boundary-Layer  Flows,"  H.L.  Reed, 
Invited  Lecture,  Tohoku  University,  Sendai,  Japan,  April  1988. 

3.  "Computational  Simulation  of  Three-Dimensional  Boundary-Layer  Flows,"  H.L.  Reed, 
Invited  Lecture,  Hokkaido  University,  Sapporo,  Japan,  April  1988. 

4.  "Stability  and  Transition  of  Compressible  Boundary  Layers,"  H.L.  Reed,  Invited  Talk, 
McDonnell  Douglas,  St.  Louis,  May  26,  1988. 

5.  "Stability  and  Transition  of  Three-Dimensional  Boundary  Layers,"  W.S.  Saric  and  H.L. 
Reed,  Invited  Paper,  AGARD  Conference  No.  438,  Fluid  Dynamics  of  Three- 
Dimensional  Turbulent  Shear  Flows  and  Transition,  Cesme,  Turkey,  Oct.  1988. 

6.  "Energy-Efficient  Aircraft,"  H.L.  Reed,  Invited  Talk,  Society  of  Women  Engineers, 
Notre  Dame,  Nov.  9,  1988. 

7.  "Three-Dimensional  Boundary-Layer  Stability,"  H.L.  Reed,  Invited  Talk,  IBM  Lecture 
Series,  Notre  Dame,  Nov.  9,  1988. 

8.  "Three-Dimensional  Boundary-Layer  Stability,"  H.L.  Reed,  Invited  Lecture,  University 
of  Western  Ontario,  London,  Canada,  Nov.  23,  1988. 

9.  "Supersonic/Hypersonic  Stability,"  H.L.  Reed,  Invited  Presentation,  NASA/Langley 
Research  Center,  Jan.  30,  1989. 

10.  "Three-Dimensional  Boundary-Layer  Stability,"  H.L.  Reed,  Invited  Seminar.  Univ 
Houston,  Mar.  2,  1989. 

11.  "Supersonic/Hypersonic  Stability  and  Transition,"  H.L.  Reed,  Invited  Presentation, 
General  Dynamics/Fort  Worth  Division,  Mar.  17,  1989. 

12.  "Supersonic/Hypersonic  Laminar/Turbulent  Transition,"  H.L.  Reed,  G.K.  Stuckert,  and 
P.  Balakumar,  Invited  Paper,  in  Developments  in  Mechanics,  Vol.  1 5,  Proceedings  of 
the  21st  Midwestern  Mechanics  Conference,  Aug.  13-16,  1989. 

13.  "Stability  of  High-Speed  Chemically  Reacting  and  Three  Dimensional  Boundary 
Layers,"  H.L.  Reed,  G.K.  Stuckert,  and  P.  Balakumar,  S  d  lUTAM  Symposiiun  on 
Laminar-Turbulent  Transition,  Sept.  1989,  Toulouse,  ed  R.  Michel  and  D.  Amal, 
Springer-Verlag,  New  York,  1991. 

14.  "Stability  of  High-Speed  Chemically  Reacting  and  Three-Dimensional  Boundary 
Layers,"  H.L.  Reed,  Invited  Seminar,  Univ  Va,  Oct.  26,  1989. 

15.  "Stability  of  High-Speed  Chemically  Reacting  and  Three-Dimensional  Boundary 
Layers,"  H.L.  Reed,  ICASE,  NASA/Langley  Research  Center,  Oct.  31,  1989. 

16.  "Stability  of  Hypersonic,  Chemically  Reacting  Boundary  Layers,"  G.K.  Stuckert  and 
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H.L.  Reed,  Bull.  Amer.  Phys.  Soc.,  Vol.  34,  No.  10,  Nov.  1989. 

17.  "Transition  in  High-Speed  Flows,"  H  L.  Reed,  Invited  Presentation,  H.L.  Reed,  General 
Dynamics,  Fort  Worth,  April  9,  1990. 

18.  "Stability  Limits  of  Supersonic  Three-Dimensional  Boundary  Layers,"  H.L.  Reed,  T. 
Haynes,  and  P.  Balakumar,  AIAA  21st  Fluid  Dynamics,  Plasmadynamics  and  Lasers 
Conference,  AIAA  Paper  90-1528,  Seattle,  June  18-20,  1990. 

19.  "Stability  of  Chemically  Reacting,  Hypersonic,  Viscous  Flows,"  G.K.  Stuckert  and  H.L. 
Reed,  AIAA  21st  Fluid  Dynamics,  Plasmadynamics  and  Lasers  Conference,  AIAA 
Paper  90-1529,  Seattle,  June  18-20, 1990. 

20.  "Stability  of  High-Speed  Boundary-Layer  Flows,"  H.L.  Reed,  T.S.  Haynes,  and  G.K. 
Stuckert,  NASA/Langley  Research  Center,  August  15, 1990. 

21.  "Crossflow  Instability  in  Supersonic  Flow,"  H.L.  Reed  and  T.S.  Haynes,  43th  Meeting 
of  the  American  Physical  Society,  Division  of  Fluid  Dynamics,  Cornell,  Bulletin 
American  Physical  Society,  Volume  35,  Number  10,  Page  2322,  November  1990. 

22.  "Spinning-Cone  Computations,"  H.L.  Reed  and  T.S.  Haynes,  Invited  Presentation,  4th 
NASP  Transition  Workshop,  NASA/Langley  Research  Center,  December  4,  1990. 

23.  "Crossflow  Calculations  for  Supersonic  Swept  Wings,"  H.L.  Reed,  Invited  Seminar, 
McDonnell-Douglas,  St.  Louis,  May  13-14,  1991. 

24.  "Stability  of  High-Speed  Flows,"  H.L.  Reed,  Eli  Reshotko's  60th  Birthday  Party, 
ICASE  and  NASA/Langley  Research  Center,  July  28, 1991. 

25.  "Stability  of  Hypersonic  Boundary-Layer  Flows  with  Chemistry,"  H.L.  Reed,  G.K. 
Stuckert,  and  T.S.  Haynes,  70th  AGARD  Fluid  Dynamics  Panel  Symposium  on 
Theoretical  and  Experimental  Methods  in  Hypersonic  Flows,  Torino,  Italy,  May  4-8, 
1992. 


Post  Doctoral  Associates 

P.  Balakumar,  "Stability  of  Three-Dimensional  Supersonic  Boundary  Layers,"  completed 
Spring  1989. 

Ph.D.  Students 

G.  Stuckert,  "Linear  Stability  Theory  of  Hypersonic,  Chemically  Reacting  Viscous  Flows," 
completed  May  1991. 

T.  Haynes,  "Effect  of  Initial  Conditions  on  Three-Dimensional  Supersonic  Boundary-Layer 
Flows,"  expected  Spring  1994. 

MS  Students 

G.  Stuckert,  "Hypersonic  Viscous  Flow  over  Two-Dimensional  and  Axisymmetric  Bodies," 
completed  Spring  1987. 

T.  Haynes,  "Stability  Limits  of  Three-Dimensional  Supersonic  Boundary-Layer  Flows," 
completed  May  1991. 

T.  (Taylor)  Trekas,  "Attachment-Line  Boundary  Layer  in  High-Speed  Flows,"  Garrett 
Fellow,  expected  Summer  1992. 

M.  Koop,  "Effect  of  Basic  State  on  Three-Dimensional  Supersonic  Boundary- Layer 
Stability,"  Garrett  Fellow,  expected  Summer  1992. 

K.  Shepherd,  "Attachment-Line  Instabilities  in  High-Speed  Flows,"  expected  Summer  1992. 

Undergraduate  Students 
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T.  Haynes,  "Stability  of  Boundary-Layer  Flows,"  completed  Fall  1988. 

T.  Taylor  (Trekas),  "Stability  of  High-Speed  Flows,"  completed  Summer  1989. 

G.  Loring,  "Secondary  Instabilities  in  Three-Dimensional  Boundary  Layers,"  completed 
Summer  1989. 

D.  Fuciarelli,  "High-Frequency  Breakdown  in  Three-Dimensional  Boundary  Layers," 
expected  Spring  1992. 

M.  Petraglia,  "Curvature  Effects  in  Three-Dimensional  Boundary-Layer  Flows,"  expected 
Spring  1993. 


The  basic  accomplishments  that  are  described  in  these  publications  are  outlined  below.  This 
work  was  invited  by  Dennis  Bushnell  to  the  4th  NASP  meeting  December  4,  1990  at 
NASA! Langley  Research  Center  as  one  of  the  few  efforts  in  3-D  high-speed  transition. 

1.  Successful  development  of  hypersonic  basic-state  Parabolized-Navier-Stokes  flow  over 
a  flat  plate  and  cone  including  nonequilibrium  chemistry 

2.  Successful  development  of  3-D  linear-stability  code  including  nonequilibrium 
chemistry  for  application  to  (I).  First  results  in  literature. 

3.  Formulate  realistic  shock  boundary  conditions  for  disturbance  state. 

4.  Comparisons  of  ideal-gas,  equilibrium  chemistry,  and  non-equilibrium  chemistry  results 
for  Mach  number  25.  First  results  in  literature. 

5.  Successful  development  of  3-D  supersonic  basic-state  boundary-layer  code  for  flow 
over  a  rotating  cone.  This  is  a  model  for  crossflow  instability. 

6.  Successful  development  of  3-D  supersonic  linear-stability  code  including  curvature  for 
application  to  (5). 

7.  Using  (6),  we  find  for  3-D  boundary  layers  that  the  second  mode  is  now  oblique  and 
cooling  is  ineffective  for  crossflow.  First  extensive  results  in  literature. 

8.  Successful  correlation  between  a  newly  formulated  crossflow  Reynolds  number  and  N 
=  9  (theoretical  transition  location)  for  the  rotating  cone.  Successful  verification  with  the 
yawed-cone  experiments  of  King  in  NASA/Langley  Mach  3.5  Quiet  Tunnel  and  of  Stetson 
at  Mach  5.9  This  is  extremely  important  to  industry. 

9.  Supersonic  attachment-line  basic-state  calculations  underway. 

10.  Successful  development  of  3-D  supersonic/hypersonic  basic-state  Parabolized-Navier- 
Stokes  code  for  flow  over  a  rotating  cone.  This  is  a  model  for  crossflow  instability  and 
includes  nonequilibrium  chemistry.  These  will  be  extended  to  real-wing  geometries. 

1 1 .  Parabolized  Stability  Equations  under  development  for  the  crossflow  instability. 
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з.  LIST  OF  SYMBOLS 

A  =Pri«(Y- l)M^2/2 

C*  Chapman-Rubesin  parameter  based  on  T* 

C%d  Chapman-Rubesin  parameter  based  on  T*j,d 

Et  total  energy  of  mixture  per  unit  volume 
h,  =1  +  Ky 

hs  =(r  +  y  cose)'" 

H  factor  in  R^fCnew)  including  compressibility 

Hcf  =  5io  /  ymaxi  crossflow  shape  factor 

i  square  root  of  - 1 

ji  molar  flux  of  species  i 

L  factor  in  Rj.f(new)  including  cooling 

m  =0  for  2-D  flow;  =1  for  axis','mmetric  flow 
M  freestream  Mach  number 

Me  local  edge  Mach  number  with  respect  to  a  reference  frame  fixed  to  surface  of  rotating 

cone 

Mr  local  Mach  number  of  basic  state  relative  to  velocity  at  generalized  inflection  points 

M<„  pre-shock  Mach  number 

N  amplification  factor 

p  pressure;  mixture  pressure 

Pr  Prandtl  number  at  the  reference  temperature  T* 

q  general  basic- state  quantity 

qjj.qy.q^heat  flux 

qo  general  disturbance-amplitude  distribution 

qi  general  disturbance  quantity 

r  distance  from  body  axis  to  surface 

r  =  T^^,o  /  Tjo 

R  Reynolds  number  based  on  local  edge  conditions  and  reference  boundary-layer  thickness 
^cf  =  Wn,„5io/V;  traditional  crossflow  Reynolds  number 
Rcf<new)  new  crossflow  Reynolds  number,  Eq.  (1) 

ReL  Reynolds  number  based  on  a  reference  length 
Ri  R  at  transition  according  to  linear  stability  theory 
Ro  R  at  Xq 

S,  Sjj  stress  tensor 

t  time 

T  temperature 

Tjd  locad  adiabatic  wall  temperature 

Tao  adiabatic  wall  temperature  at  upstream  point  Xq 

Tj  edge  temperature 

T^o  specified  wall  temperature  at  upstream  point  Xq 

local  wall  temperature 
T*  reference  temperature 

T*ad  reference  temperature  for  adiabatic  wall 

и, v,w  streamwise,  normal-to-wall,  spanwise  velocity 

Ug  local  inviscid  flow  speed;  for  rotating  cone,  reference  frame  is  fixed  to  the  surface 

maximum  crossflow  velocity,  that  is,  velocity  in  direction  perpendicular  to  local  inviscid 
flow;  for  rotating  cone,  reference  frame  is  fixed  to  surface 
x,y,z  streamwise,  normal-to-wall,  spanwise  coordinates 

Xq  initial  streamwise  position  for  calculation  of  N;  or  upstream  position  for  marching 

Ymax  height  of  W„,„  above  wall 

yi  height  of  relative  sonic  point  above  wall 

a  streamwise  component  of  wavenumber 

-tti  negative  of  imaginary  part  of  a;  spatial  growth  rate 
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P  spanwise  component  of  wavenumber 

Pi  imaginary  part  of  P 

S  boundary-layer  thickness 

5ad  boundary-layer  thickness  for  adiabatic  wall 
5cooi  boundary-layer  thickness  for  cooled  wall 

Sincomp  boundary-layer  thickness  for  incompressible  constant-temperature  flow 

5r  =  (ite^PeUe)'^;  reference  boundary-layer  thickness 

5io  height  above  y^^^j  where  crossflow  is  10% 

e  angle  between  body  surface  and  body  axis 

Y  edge  ratio  of  specific  heats 

Yoo  pre-shock  ratio  of  frozen  specific  heats 

r|  either  a  similarity  or  computational  variable  in  normal-to-wall  direction 

K  =  -de  /  dx 

Pe  edge  dynamic  viscosity 

Vg  edge  kinematic  viscosity 

9  cone  angle 

p  mixture  mass  density 

pe  edge  density 

Oj  ratio  of  species  i  concentration  to  mixture  mass  density 

T  time  variable  in  computational  domain 

(0  frequency  of  disturbance 

Q  cone  rotational  speed 

streamwise,  normal-to-wall,  spanwise  computational  coordinates 

4.  MOTIVATION 

The  skin-friction  drag  and  heat-transfer  rates  of  hypersonic  and  supersonic  vehicles  can 
be  reduced  by  delaying  transition  from  laminar  to  turbulent  flow.  But  the  physical 
characteristics  of  transition  at  these  speeds  are  not  well  known.  Transition  is  known,  however,  to 
be  highly  dependent  on  the  details  of  the  flowfield. 

Hypersonic  and  supersonic  flows  are  further  complicated  over  incompressible  flows  for 
some  of  the  following  reasons.  1)  At  these  speeds,  the  gas  often  cannot  be  modeled  as  perfect 
because  the  molecular  species  begin  to  dissociate  due  to  aerodynamic  heating.  In  fact, 
sometimes  ^here  are  not  enough  intermolecular  collisions  to  support  local  chemical  equilibrium 
and  a  nonequilibrium-chemistry  model  must  be  used.  2)  The  bow  shock  is  very  close  to  the 
edge  of  the  boundary  layer  and  must  be  included  in  the  calculations.  3)  Because  of  the  speeds, 
the  wings  are  highly  swept  and  the  boundary  layers  are  highly  3-D.  These  effects  must  be 
included  in  any  studies  of  transition. 

Linear  stability  theory  proves  useful  in  determining  the  important  effects  and  their  trends. 
Here,  the  resultant  growtit  or  decay  of  small  disturbances  in  the  boundary  layer  which  lead  to 
transition  is  the  measure  of  transition  enhancement  or  delay.  This  paper  reports  the  results  of 
linear  stability  theory  applied  to  chemically  reacting  flows  with  finite  shock  stand-off  distance  in 
Section  5  and  3-D  flows  in  Section  6. 

5.  CHEMICALLY  REACTINXJ  TWO-DIMENSIONAL  BOUNDARY  LAYERS 
5.1  Numerical  Approach 

Qualitative  as  well  as  quantitative  differences  in  the  linear  stability  of  the  shock  layer  on 
a  sharp  cone  can  be  observed  when  the  effects  of  finite-rate  chemical  reactions  are  included  in 
the  analysis.  Here  a  five-component  model  for  dissociated  air  is  used:  O2,  N2,  NO,  N,  and  O. 
Only  the  effects  of  dissociation  are  considered;  those  of  ionization  are  not.  The  mixture  is  also 
assumed  to  be  one  of  ideal  gases  in  thermal  equilibrium. 

The  equations  governing  the  flow  of  this  reacting-gas  mixture  are  originally  expres.sed  in 
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the  body-intrinsic  coordinate  system  (i.e.,  Anderson  et  al.,  1984).  However,  since  the  boundary 
layer  occupies  a  substantial  fraction  of  the  shock  layer  at  hypersonic  speeds,  these  equations  are 
transformed  into  a  coordinate  system  where  the  basic-state  bow  shock  is  a  boundary  of  the 
computational  domain; 

3(y,iQ)/3t  +  a(yT,E)/a^  +  af-y^Q  -  y^E  +  F  -y;G)/an  +  a(y^G)/ai;  +  y^H  -  y^R  =  0 

where 

Q  =  hih3  {p,  pu,  pv,  pw,  Ej,  pOi) 

E  =  Ej  -  Ey 

Ej  =  h3  {pu,  pu24-p/Y„M„2,  puv,  puw,  (E,+p)u,  pCjU) 

Ey  =  (h3/ReL)  {0,  S„,  S^,  YooM„2  (uS„  +  vSy,  -i-  wS„)  -  q^,  -ji,) 

F  =  Fj-Fy 

Fj  =  hih3  (pv,  puv,  pv2+p/YJVl„2,  pyw,  (Ei-i-p)v,  pOiV} 

Fy  =  (hlh3/ReL)  {0,  S,y,  Syy,  S,y,  YocM„2  (uS^y  +  VSyy  -1-  WS^)  -  qy,  -jjy} 

G  =  Gf  -  Gy 

Gi  =  h]  (pw,  puw,  pvw,  pw2-i-p/Y„M„2  (E^+p)w,  po^w) 

Gy  =  (h,/ReO  (0,  Sy„  y»M„2  (uS„  +  vSy^  +  wS  J  -  q^,  -jjj 

H  =  Hi  -  Hy 

Hi  =  (0,  Kh3puv  -  mh,  (pw2  +  sine,  -Khj  (pu2  +  -  mh]  (pw2  +  p/Yo„M„2)  cose, 

mh|pw  (usine  +  vcose),  0, 0} 

Hy  =  (0,  Kh3Sxy  -  mhiS^sine,  -KhjS^  -  mhiSj,^cose,  mhi  (S;^2sine  +  Sy^cose),  0, 0) 

R  =  hih3  (0,  0,  0,  0,  0,  Molar  rate  of  production  of  species  /} 

(Note;  E,  F,  G,  H,  Q,  R  are  not  in  the  List  of  Symbols.)  The  viscosity  used  to  determine  the 
viscous  stress  tensor  S  is  computed  using  the  mixture  laws  of  Brokaw  (1958).  The  translational 
thermal  conductivity  is  computed  similarly,  whereas  the  internal  thermal  conductivity  is 
determined  as  described  in  Hirschfelder  (1957).  The  molar  fluxes  jj  are  computed  using  the 
multicomponent  diffusion  model  described  by  Cuniss  &  Hirschfelder  (1949),  but  only  diffusion 
due  to  concentration  gradients  is  included  -  diffusion  due  to  pressure  gradients,  temperature 
gradients,  and  body  forces  is  not.  Finally,  the  law  of  mass  action  is  used  to  compute  the  molar 
rate  of  production  of  each  species  assuming  that  they  participate  in  the  following  elementary 
reactions; 

02  +  M<->0  +  0  +  M 
N2  +  M<->N  +  N  +  M 
N2  +  O  <->  NO  +  N 
O  +  NO  N  +  O2 

where  M  is  a  collision  partner  (any  of  the  species  present  in  the  mixture)  which  transfers  energy 
in  a  reaction. 

The  thermophysical  data  needed  for  the  analysis  have  been  taken  from  a  variety  of 
sources.  Collision  cross  section  data  for  the  transport  properties  can  be  found  in  Biolsi  &  Biolsi 
(1983),  Biolsi  (1988),  Capitelli  &  Devoto  (1973),  Capitelli  &  Ficocelli  (1972),  Cubley  &  Mason 
(1975),  Levin  et  al.  (1987,  1988),  Monchick  (1959),  and  Yun  &  Mason  (1962).  Thermodynamic 
data  have  been  taken  from  Blottner  et  al.  (1971)  and  Jaffe  (1987).  Reaction-rate  data  can  be 
found  in  Camac  &  Vaughan  (1961),  Wray  (1962),  Thielen  &  Roth  (1986),  Monat  et  al.  (1978), 
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and  Hanson  &  Salimian  (1984).  Reverse  reaction-rate  constant.',  are  computed  using  the  law  of 
detailed  balance  express  them  in  terms  of  the  forward-rate  constants  and  the  equilibrium 
constants.  For  a  complete  description  of  the  constitutive  equations  as  well  as  detailed  references 
for  the  thermophysic^  data,  refer  to  Stuckert  (1991). 

5.1.1  Basic  State 

The  basic  state  is  computed  using  the  Parabolized  Navier-Stokes  equations,  which  are 
obtained  from  those  given  above  by  neglecting  viscous  terms  involving  derivatives  tangent  to  the 
body.  (To  simplify  the  analysis,  this  approximation  is  made  in  the  body-intrinsic  coordinate 
system  before  the  equations  are  transformed  into  the  shock-fit  coordinate  system.)  The  Beam- 
Warming  algorithm  (Beam  &  Warming,  1978)  is  used  in  conjunction  with  a  shock-fitting 
scheme  and  the  Vigneron  technique  (to  suppress  departure  solutions)  to  solve  these  equations. 
This  is  an  efficient  implicit  algorithm  used  to  march  the  solution  downstream  from  initial 
conditions  generated  here  by  assuming  that  the  flow  is  approximately  conical  near  the  apex  of  a 
sharply  pointed  body.  A  stretching  function  is  used  to  cluster  grid  points  near  the  edge  of  the 
boundary  layer  -  a  location  where  the  basic  state  varies  rapidly  and  which  also  lies  in  the  vicinity 
of  the  critical  point  of  the  disturbance  state  at  hypersonic  Mach  numbers.  No  artificial 
dissipation  is  required  for  the  simple  geometry  investigated  here  because  there  are  no 
discontinuities  embedded  in  the  shock  layer. 

5.1.2  Linear  Stability  Theory 

The  linear  stability  analysis  for  each  of  these  cases  is  performed  by  linearizing  the 
complete  Navier-Stokes  equations  about  the  basic  state.  The  same  basic-state  shock-fit 
coordinate  system  is  used  (which  is  almost  conical)  and  linearized  inviscid  shock  jump 
conditions  are  imposed  at  the  basic-state  shock  location.  For  a  cone,  though,  k  s  0,  and  so  hj  s 
1.  Normal  modes  are  assumed  for  the  disturbance  state  variables  qj: 

qi  =  qo(^)  exp{i(a^  -h  PC  -  on)} 

and  the  basic-state  flow  is  assumed  to  satisfy  dq/d^  =  0.  (Note  that  the  basic-state  normal- 
velocity  component  is  not  set  to  zero;  it  is  not  zero  in  a  conical  flow.)  The  disturbances  in 
velocity  and  temperature  as  well  as  species  mass  fluxes  (i.e.,  noncatalytic  wall)  are  zero  at  the 
surface.  (For  equilibrium  air  analyses,  the  disturbances  are  assumed  to  be  in  chemical 
equilibrium  at  the  surface.)  These  homogeneous  equations  and  boundary  conditions  represent  an 
eigenvalue  problem  for  the  eigenvector  qofq).  When  a  and  p  are  specified,  the  eigenvalue  is  to 
(temporal  problem),  and  when  P  and  co  are  specified,  the  eigenvalue  is  a  (spatial  problem).  In 
either  case,  the  governing  equations  are  first  discretized  using  second-order  accurate  central 
differences.  The  temporal  problem  is  easy  to  analyze  because  it  is  linear  in  co  and  hence  can  be 
solved  using  a  globally  convergent  scheme  such  as  the  QZ  algorithm  (e.g.,  Golub  &  VanLoan, 
1985).  Inverse  iteration  and  the  Newton-Raphson  method  can  then  be  used  to  compute  the 
eigenvector  and  solve  the  spatial  eigenvalue  problem. 

5.2  Results 

The  particular  case  investigated  is  that  of  a  10°  half-angle  cone  flying  at  81(X)  m/s  and 
zero  angle  of  attack  at  an  altitude  of  60.96  km  where  the  ambient  temperature  and  pressure  are 
252. 6K  and  2.008X10'*  atmospheres,  respectively.  The  Reynolds  number  based  upon  a 
reference  length  of  Im  and  these  freestream  conditions  is  139,900.  The  atmosphere  is 
considered  to  be  21%  O2  and  79%  N2  by  vo' ame,  and  the  noncatalytic  surface  of  the  cone  is  kept 
at  a  constant  temperature  of  1200K.  This  is  essentially  the  case  considered  by  Prabhu  et  al. 
(1988)  except  that  they  assumed  a  different  freestream  composition  of  air.  Also,  there  is  no 
crossflow. 

Due  to  the  lack  of  space,  only  a  couple  of  basic-state  results  pertinent  to  the  stability 
analysis  are  presented  here.  It  should  be  mentioned,  though,  that  the  boundary-layer  edge  Mach 
number  is  approximately  10  for  the  given  pre-shock  conditions.  It  is  slightly  higher  for  the 
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perfect-gas  analysis,  and  slightly  lower  for  the  equilibrium-  and  nonequilibrium-air  analyses 
because  the  hotter  perfect-gas  boundary  layer  is  more  rarefied  and  hence  thicker,  displacing  and 
compressing  the  external  flow  the  greatest.  Similarly,  the  perfect-gas  edge  temperature  is 
approximately  6.5  times  the  freestream  temperature,  whereas  the  equilibrium-  and 
nonequilibrium-air  edge  temperatures  are  closer  to  6  times  the  freestream  temperature. 


One  feature  of  the  boundary  layer  important  to  the  linear  stability  is  the  size  of  the  region 
where  the  disturbance  state  is  supersonic  relative  to  the  local  basic  state.  A  measure  of  this  is 

yi 

(M,2-l)i«dy 


J  0 


where  is  the  local  Mach  number  of  the  basic  state  relative  to  the  velocity  at  the  generalized 
inflection  points  (those  points  where  3(p9u/3y)/8y  =  0).  This  integral  is  shown  in  Figure  1 
(where  x  is  measured  in  meters).  Corresponding  results  from  a  perfect-gas  boundary-layer 
similarity  solution  are  also  shown  for  comparison.  One  can  see  from  this  figure  that  the  perfect- 
gas  similarity  solution  and  PNS  solution  yield  very  comparable  results  and  that  in  both  cases  the 
values  of  the  integrals  are  substantially  less  than  the  corresponding  values  for  the  equilibrium- 
and  nonequilibrium-air  boundary  layers.  As  Mack  (1969)  has  shown,  this  integral  is  closely 
related  to  the  spacing  between  his  higher  inviscid  modes;  as  an  approximation,  it  is  inversely 
proportional  to  it.  The  shift  to  lower  frequencies  of  the  most  unstable  second-mode  disturbance 
when  the  effects  of  equilibrium-air  chemistry  are  included  has  been  seen  by  both  Malik  (1989) 
and  Stucken  (1991).  See  Figure  2.  The  same  shift  in  frequency  is  observed  in  the 
nonequilibrium  case  for  precisely  the  same  reason;  the  local  relative  Mach  ;.umber  is  in  general 
higher  for  the  reacting-gas  cases  because  the  local  velocity  is  almost  the  same,  but  the  local 
speed  of  sound  is  lower  in  the  colder  boundary  layers.  These  reacting-gas  boundary  layers  are 
colder  because  some  of  the  heat  generated  by  viscous  dissipation  is  used  to  dissociate  the 
molecular  oxygen  and  nitrogen. 

The  equilibrium-air  boundary  layer  also  exhibits  a  feature  which  is  not  seen  in  the  perfect 
gas  or  the  nonequilibrium-air  boundary  layer  (at  least  for  the  present  flight  conditions).  In 
particular,  an  equilibrium-air  mode  exists  which  is  supersonic  relative  to  the  inviscid  region  of 
the  shock  layer.  The  eigenfunction  for  this  disturbance  is  shown  in  Figure  3.  The  oscillatory 
behavior  of  the  magnitude  of  the  disturbance  is  due  in  part  to  the  fact  that  the  character  of  ihe 
disturbance-state  equations  changes  when  the  Mach  number  of  the  disturbance  is  supersonic 
relative  to  the  boundary- layer  edge.  It  is  also  due  to  the  fact  that  the  shock  is  a  finite  distance 
av/ay  from  the  surface.  Since  the  shock  stand-off  distance  is  finite,  solutions  to  the  disturbance- 
state  equations  in  the  inviscid  region  of  the  shock  layer  are  possible  v»hich  do  not  decay  .,s  y 
becomes  large.  This  is  discussed  in  greater  detail  by  Stuckert  &  Reed  (1992). 


6.  THREE  -  DIMENSIONAL  HIGH  -  SPEED  BOUNDARY  LAYERS 
6.1  Use  of  Simple  Geometries  to  Model  Three-Dimensional  Boundary  Layers 

At  hypersonic  and  supersonic  speeds,  the  wings  of  the  vehicles  arc  highly  swept  and  the 
flow  is  highly  3-D.  When  a  boundary-layer  flow  is  fully  3-D,  the  stability  and  transition 
behav'or  is  quite  different  from  2-D  flows.  Of  particular  interest  are  the  stability  characteristics 
of  these  3-D  flows  where  inviscid  criteria  may  produce  a  stronger  instability  (crossflow 
instability)  than  the  usual  Tollmien-Schlichting  waves  (T-S  waves).  For  a  recent  review  of  this 
subject,  see  Reed  &  Saric  (1989). 

Examples  of  other  3-D  flows  of  practical  interest  include  rotating  cones,  yawed  cylinders 
and  cones,  comers,  inlets,  and  rotating  disks.  It  appears  that  all  of  these  flows  (including  swept 
w'ngs)  exhibit  the  same  rich  variety  of  stability  behavior  that  is  generic  to  3-D  boundary  layers 
A  consistent  characteristic  of  the  instabilities  is  the  presence  of  streamwise  vorticity  within  the 
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shear  layer. 

In  fact,  in  incompressible  flow,  for  example,  the  experimental  and  analytical  results  of 
Kobayashi  et  al.  (1983)  on  rotating  cones  have  demonstrate  that  one  can  even  model  a  swept 
wing  or  other  more  complicated  geometry  exhibiting  crossflow  instability  by  a  simpler  geometry 
and  study  the  stability  characteristics.  Results  obtained  can  thus  be  directly  applied  to  other  3-D 
boundary  layers.  From  simpler  geometries  such  as  the  rotating  cone,  much  of  what  is  known 
about  the  flow  over  the  swept  wing  has  been  learned. 

6.2  Numerical  Approach 

To  evaluate  parameters  quantifying  stability  characteristics,  the  linear  stability  of  the 
flow  over  a  rotating  cone  at  zero  incidence  is  examined.  As  mentioned  above,  this  is  a  simple 
geometry  that  is  often  used  very  successfully  in  incompressible  flow  as  a  model  for  a  swept 
wing.  At  very  high  speeds,  where  even  the  basic-state  calculations  are  a  problem,  the  simple 
geometry  of  the  rotating  cone  becomes  a  suitable  and  valuable  model  to  study  the  crossflow 
instability  (3alakumar  &  Reed,  1991). 

The  geometry  considered  is  that  of  a  sharp  cone  located  in  a  uniform  high-speed  stream 
at  zero  angle  of  attack  and  rotating  about  its  axis  (Figure  4).  The  governing  boundary-layer 
equations  for  a  compressible  ideal-gas  flow  are  solved  in  a  body-oriented  coordinate  system.  By 
varying  the  freestream  Mach  number,  rotational  speed,  cone  angle,  and  position  on  the  cone,  a 
wide  parameter  range  of  possible  3-D  boundary  layers  can  be  studied  and  the  various  non- 
dimensional  parameters  associated  with  the  crossflow  profile  are  easily  fixed  (R,  R(.f,  H^.^,  .  .  .). 
The  details  of  the  basic-state  formulation  and  linear  stability  analysis  are  available  in  Balakumar 
&  Reed  (1991). 

6.3  Linear  Stability  Theory  Results 

Calculations  were  completed  at  different  edge  Mach  numbers  M  =  1.5,  3.0,  5.0  and  8.0 
and  Reynolds  numbers  R  =  600,  1000,  2000,  3000  for  a  cone  of  half  angle  15°  at  a  non- 
dimensional  rotational  speed  of  0.375  by  Balakumar  &  Reed  (1991).  To  summarize,  they  find  at 
all  Mach  numbers  that  the  maximum  amplified  crossflow  is  nonstationary  with  the  frequency 
increasing  with  Mach  number  from  approximately  5  kHz  at  M=1.5.  There  exist  neutral  and 
unstable  stationary  waves  but  the  amplification  rate  is  very  small  compared  to  the  positive 
frequency. 

At  M  =  1.5,  the  direction  of  the  most  amplified  crossflow  is  55°  at  R  =  600  increasing  to 
65°  at  R  =  3000.  As  R  increases  further  the  angle  approaches  that  of  the  crossflow  direction  of 
73.8°.  A  similar  observation  occurs  for  other  Mach  numbers. 

At  M  =  5.0,  for  which  the  2-D  second  mode  has  been  observed  to  clearly  dominate  the 
first  mode  in  the  flat-plate  boundary  layer,  they  find  a  competition  between  the  crossflow  and 
second  mode.  At  R  =  600,  the  maximum  amplification  rate  of  the  crossflow  is  higher  than  the 
maximum  amplification  rate  of  the  2-D  second  mode;  see  Figure  5  (from  Balakumar  &  Reed, 
1991).  As  Reynolds  number  increases,  even  though  the  amplification  rate  of  the  crossflow 
increases,  it  is  eventually  less  than  the  maximum  amplification  of  the  second  mode.  As  Mach 
number  is  increased,  the  second  mode  dominates  earlier. 

Further  results  show  that  the  most  unstable  second  mode  in  a  3-D  boundary  layer  is 
actually  oblique  whereas  the  second  mode  in  a  2-D  boundary  layer  is  2-D;  this  is  also  shown  in 
Figure  5  where  the  2-D  direction  is  -16.2°.  The  most  unstable  second  mode  is  inclined  at  5°  to 
the  invi.scid  flow  direction  at  M  =  5  and  is  inclined  at  9°  at  M  =  8.  However,  the  difference 
between  the  amplification  rates  in  the  most  amplified  direction  and  in  the  inviscid  flow  direction 
is  small,  on  the  order  of  2%. 

Cooling  the  wall  does  not  affect  the  crossflow  instability.  For  M  =  5,  Figure  6  (from 
Balakumar  &  Reed,  1991)  shows  the  variation  of  the  critical  Reynolds  number  (at  which  an 
instability  first  appears)  versus  crossflow  Reynolds  number  (achieved  by  varying  rotational 
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speed,  a  larger  value  corresponds  to  stronger  three-dimensionality)  for  different  wall  conditions  r 
=  T^c/Tjo  =  1-0,  0.8,  and  0.5.  As  expected,  the  critical  Reynolds  number  of  the  first  mode 
increases  and  the  critical  Reynolds  number  of  the  second  mode  decreases  with  wall  cooling  and 
confirms  the  result  that  wall  cooling  stabilizes  the  first  mode  and  destabilizes  the  second  mode. 
However,  for  the  first  mode,  it  is  seen  that  the  critical  Reynolds  number  decreases  with 
increasing  crossflow  Reynolds  number  and  at  large  crossflow  Reynolds  number,  the  critical 
Reynolds  numbers  all  approach  the  same  value  regardless  of  the  status  of  the  boundary  layer. 
This  is  observed  for  all  Mach  numbers  and  Reynolds  numbers  tested.  Similar  results  have  been 
seen  by  Lekoudis  (1980),  Mack  (1980),  and  Bushnell  &  Malik  (1987). 

6.4  Transition  Prediction  and  Correlation  Parameters 

The  state-of-the-art  transition-prediction  method  still  involves  linear  stability  theory 
coupled  with  an  e^  transition-prediction  scheme  (Mack,  1984;  Poll,  1984)  and  is  applied  at  all 
speeds  (Bushnell  et  al.,  1989).  N  is  the  result  of  the  integration  of  the  linear  growth  rate  from  the 
first  neutral-stability  point  to  a  location  somewhere  downstream  on  the  body.  Thus  e^'  is  the 
ratio  of  the  amplitudes  at  the  two  points  and  the  method  correlates  the  transition  Reynolds 
number  with  N. 

Aside  from  the  e^  method  along  with  various  modifications  (e.g.,  Cebeci  et  al.,  1988), 
several  past  investigators  have  identified  non-dimensional  parameters  quantifying  the 
characteristics  of  the  crossflow  instability  for  boundary-layer  profiles  and,  in  some  cases, 
attempted  to  correlate  these  numbers  with  transition  location.  Some  examples  follow. 

In  subsonic  flow,  Owen  &  Randall  (1952)  introduced  a  crossflow  Reynolds  number  R^f  = 
W^„5i(/Ve  as  the  governing  parameter  and  suggested  that  transition  occurs  when  the  crossflow 
Reynolds  number  becomes  equal  to  175.  Pfenninger  (1977)  used  crossflow  shape  factor  H^f  = 
ymax/^io  crossflow  Reynolds  number  in  the  design  of  supercritical  airfoils.  Dagenhart  (1981) 
then  considered  stationary  crossflow  vortices  and,  instead  of  solving  the  linear-stability 
equations  each  time,  he  used  a  table  lookup  of  growth  rates  based  on  the  profile  characteristics: 
crossflow  shape  factor  and  crossflow  Reynolds  number.  He  then  integrated  these  interpolated 
values  to  obtain  N-factors.  He  reported  that  his  code  MARIA  adequately  reproduces  the  stability 
results  of  the  more  complicated  stability  codes  using  the  same  physical  disturbance  model  while 
using  less  than  two  percent  of  the  computer  time. 

In  supersonic  flow.  Chapman  (1961)  and  Pate  (1978)  made  similar  conclusions  that 
crossflow  Reynolds  number  correlates  well  with  transition  location.  On  a  yawed  cone.  King 
(1991)  found  that  there  was  no  correlation  with  the  traditional  definition  of  crossflow  Reynolds 
number.  However,  when  he  reformulated  this  parameter  to  include  both  compressibility  and 
ya wed-cone  geometry  effects,  he  found  a  good  correlation  for  both  his  and  Stetson's  (1982)  data. 

6.4.1  Numerical  Approach 

With  the  current  interest  in  high-speed  flight,  there  is  also  a  keen  desire  to  determine 
correlating  parameters,  based  purely  on  basic-state  profiles,  that  can  be  easily  incorporated  into 
existing  basic-state  codes  and  will  accurately  predict  transition  location  (or  trends)  for 
crossflow-dominated  problems.  To  evaluate  appropriate  parameters  quantifying  stability 
characteristics,  the  linear  stability  of  the  How  over  a  rotating  cone  at  zero  incidence  is  again 
examined.  Results  are  then  applied  to  available  experimental  data  on  other  geometries. 

Locally,  for  a  given  frequency,  the  maximum  growth  rate  -ttj  is  found  and  Pi=0.  For  each 
frequency,  then,  the  amplification  factor  N  is  determined  by  integrating  the  growth  rate  in  the 
streamwise  direction  x  from  the  Branch  1  location  Xq.  With  the  Reynolds  numbers  at  Xq  and  x 
being  R^  and  R,  respectively, 


Reed  &  Saric:  Hypersonic  Stability  (AFOSR) 


page  14 


N  =  -2 


R 


ttidR. 


Ro 


To  evaluate  transition  location,  all  possible  frequencies  are  sampled  to  determine  where  each 
individual  disturbance  achieves  N  =  9.  The  most  upstream  of  this  locus  of  streamwise  locations 
is  deemed  the  estimated  transition  location.  The  values  of  other  non-dimensional  parameters 
characterizing  crossflow  may  then  be  evaluated  against  the  results  of  N  =  9. 

6.4.2  Compressibility  Effects 

White  (1974)  points  out  that  to  estimate  the  boundary-layer  thickness  for  a  flat  plate,  one 
should  consider  the  similarity  variable  so  that 

^Ti(5) 


5  R  /  X  = 


(Tn'e)dTl 


0 


To  extend  this  to  a  more  general  situation,  then,  a  3-D  compressible  boundary-layer,  5,  should 
thicken  with  respect  to  the  corresponding  incompressible  layer,  according  to 


6/5 


incomp 


ti(6) 

(TyTe)  dtl  /  Tl 

0 


Looking  at  the  right  hand  side  of  the  last  equation,  5inj.on,p  is  a  constant-temperature 
incompressible  value.  The  estimated  thickness  then  of  a  cooled-waJl  boundary  layer,  b^ooi,  to 
that  of  an  adiabatic  wall,  5jd,  is  (White  1974,  Reed  &  Haynes  1992). 

5cooi  /  Sad  =  (C*  /  (3-279  1.721  [T,,  /  TJ  [1  +  A] 

+  0.664  A)  /  (5  +  2.385  A) 


where 


A  =  Pri/2  (y-  1)  M,2/2 

C*  =  (T*  /  TJi-  (1  -Hi  10.4  /  Te)  /  ([T*  /  TeJ  +  1 10.4  /  T,) 

C-ad  =  (T*ad  /  Te)'^  ( 1  +  1 10.4  /  T^)  /  {[T,^  /TJ+l  10.4  /  Te) 

T* /T,  =  0.5  +  0.5  T,, /Te  +  A  /  6 
T‘ad  /  Te  =  0.5  +  0.5  (1  +  A)  -t-  A  /  6 

and  the  Chapman-Rubesin  parameter  is  approximated  as  a  constant  C*  across  the  boundary  layer. 
Because  the  incompressible  crossflow  Reynolds  number  correlates  reasonably  well  for 
incompressible  flows,  a  new  general  definition  for  R^f,  then,  is 

ReKnew)  =  H  L  Ref  =  H  L  W,^,  5,o  /  (1) 


where 


ri(5,o) 

(meldT) 


J  0 

L  =  (C*  /  C*3d)'^  (3.279  +  1.721  [T,  /  TJ  (1  +  A] 
+  0.664  A)  /  (5  +  2.385  A) 


H=ti(8,o)/ 
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and  all  temperatures  are  in  K  and  ti(5io)  is  the  value  of  "n  at  5,0-  The  quantity  Sjq  has  been  scaled 
back  to  an  incompressible  value  with  the  inclusion  of  the  new  factor  H.  To  compensate  for  a 
cooled  wall,  then,  the  factor  L  is  used.  Note  that,  because  one  is  taking  the  ratio  of  two  normal- 
to-the-wall  lengths,  the  new  factor  H  is  easily  computed  no  matter  what  the  scaling  is  for  the 
normal-to-the-wall  coordinate.  Note  also  that  R^ffnew)  reduces  to  Rj,f  for  incompressible, 
adiabatic- wall  flows. 

6.4.3  Results 

Using  the  traditional  formulations  for  crossflow  Reynolds  number  and  shape  factor 

=  Wnjajj  5io  /  V  (2) 

^cf  ~  Ymax  !  ^10’ 

Figure  7  shows  the  attempt  at  transition  correlations  for  an  upstream  Mach  number  of  M  =  3  and 
various  cone  angles  0,  rotational  speeds  Q,  and  wall  and  freestream  temperatures  T,,,  and  T^.  The 
various  flow  conditions  represented  are  documented  in  Table  1.  The  spread  in  R^j  is  on  the  order 
of  200%  and  is  therefore  not  useful  for  design.  Figure  8  shows  the  same  data  plotted  considering 
transition  Reynolds  number  Rt  versus  rotational  speed  Q  (with  increasing  Q  implying  increasing 
3-D  effects).  Cooling  is  only  slightly  stabilizing  (see  also  Lekoudis  1980,  Mack  1980,  Bushnell 
&  Malik  1987,  and  Balakumar  &  Reed  1991)  and  an  increase  of  stagnation  temperature  has  an 
even  smaller  stabilizing  effect.  Even  so,  these  temperature  effects  produce  large  changes  in  R^f. 


Table  1  Rotating-cone  configurations  for  M  =  3 


Set 

Fixed 

Varied  Parameters 

Parameters 

and  Range 

1 

a=0.375 
T,=70K 
adiabatic  wall 

0  from  10°  to  35° 

II 

Q=0.25 
Te=70K 
adiabatic  wall 

0  from  10°  to  20° 

III 

0=15“ 
Te=70K 
adiabatic  wall 

Q  from  0.1  to  0.8 

IV 

0=15“ 
T,=260K 
adiabatic  wall 

Q  from  0.1  to  0.8 

V 

(2=0,375 

0  from  10°  to  35° 

* 

T,=260K 
adiabatic  wall 

VI 

0=15“ 
T,=70K 
cooled  wall 

Q  from  O.I  to  0.8 

VII 

0=15“ 
T,=260K 
cooled  wall 

Q  from  0.4  to  0.8 

Table  2  shows  comparative  results  for  upstream  freestream  Mach  numbers  of  M  =  0.01, 
3,  and  6  for  different  ffeestream-temperature  (TJ  and  surface-temperature  (TJTJ  conditions 
and  a  wide  variety  of  rotating-cone  geometries.  The  Mach-3  (M  =  3)  conditions  correspond  to 
those  in  Figure  7  and  Table  1.  The  results  of  the  new  definition  for  crossflow  Reynolds  number 
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with  compressibility  effects  included,  Rcf(new)  [Eq.  (1)],  are  contrasted  with  those  obtained  from 
the  traditional  definition,  Rj.f  [Eq.  (2)].  Also  included  is  the  maximum  crossflow  velocity  in 
%Ue;  the  significance  of  this  should  become  apparent  in  the  subsequent  discussion. 

Figure  9  shows  all  the  rotating-cone  data  from  Table  2  for  M  =  0  and  3  plotted.  There 
appears  to  be  a  relationship  between  maximum  crossflow  velocity  and  the  new  crossflow 
Reynolds  number,  Eq.  (1),  proposed  above.  Also,  the  new  compressible  values  are  consistent 
with  the  incompressible  value. 


Table  2 


M 

Me 

Te(K) 

T  /T 

Ref 

Rcf(new)W„,„, 

.01 

.01 

300 

1 

165 

165 

5.9 

3 

3.1 

70 

2.6 

241 

119 

3.2 

3 

3.2 

70 

2.7 

311 

149 

4.5 

3 

3.4 

70 

2.9 

373 

170 

5.7 

3 

3.6 

70 

3.2 

409 

175 

6.1 

3 

3.8 

70 

3.5 

428 

171 

6.1 

3 

3.1 

260 

2.5 

210 

107 

2.6 

3 

3.1 

260 

2.6 

263 

132 

3.6 

3 

3.2 

260 

2.7 

316 

154 

4.8 

3 

3.4 

260 

2.9 

366 

168 

5.9 

3 

3.7 

260 

3.2 

388 

166 

6.1 

3 

3.9 

260 

3.5 

400 

161 

6.0 

3 

3.2 

70 

1.5 

339 

177 

4.6 

3 

3.3 

70 

1.5 

354 

179 

5.5 

3 

3.6 

70 

1.5 

354 

168 

6.1 

3 

3.8 

70 

1.5 

344 

155 

6.1 

3 

4.1 

70 

1.5 

332 

140 

5.8 

3 

4.2 

260 

1.5 

448 

201 

5.6 

3 

3.8 

260 

1.5 

323 

156 

6.0 

3 

4.1 

260 

1.5 

310 

142 

5.7 

6 

9.2 

70 

14.8 

2483 

255 

5.0 

6 

9.2 

260 

17.1 

2176 

217 

5.0 

6 

7.4 

70 

9.7 

2101 

309 

6.1 

6 

7.3 

70 

4.4 

1442 

247 

6.1 

Because  linear  stability  theory  with  N  =  9  applied  to  a  rotating  cone  was  used  to  find  this 
trend,  it  is  important  to  verify  it  against  experimental  data.  The  most  reliable,  available,  high¬ 
speed  3-D  transition  data  to  these  authors’  knowledge  is  that  of  King  (1991)  on  the  yawed  cone 
in  the  Mach  3.5  Quiet  Tunnel  at  NASA/Langley.  The  King  experiment  was  on  a  5°  half-angle 
cone  yawed  at  0.6°,  2°,  and  4°.  The  data  of  transition  locations  for  various  freestream  conditions, 
both  quiet  and  noisy,  are  documented  in  King  (1991). 

Later,  King  provided  the  present  authors  with  the  computational  mean-flow  profiles  he 
used  in  the  analysis  of  his  experimental  results.  Applying  the  traditional  definition  of  crossflow 
Reynolds  number,  Eq.  (2),  gives  values  ranging  from  80  to  640  for  quiet  conditions  and  60  to 
560  for  noisy  conditions.  King  found  a  correlation  for  his  data  when  both  a  geometry  and 
compressibility  correction  were  applied.  The  new  parameter  proposed  in  the  present  paper 
contains  no  geometry  factor  explicitly.  Applying  the  new  crossflow  Reynolds  number,  Eq.  (1), 
to  these  same  profiles  results  in  Figure  10. 

King  (1991)  also  considered  the  M  =  5.9  experiments  of  Stetson  (1982)  on  a  8°  half-angle 
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cone  yawed  at  1°,  2°,  and  4°  and  again  provided  the  present  authors  with  the  mean-flow  profiles 
he  used  in  the  evaluation  of  this  experiment.  Considering  what  Stetson  terms  as  the  beginning  of 
transition,  applying  the  traditional  definition  of  crossflow  Reynolds  number,  Eq.  (2),  gives 
values  ranging  from  140  to  780.  The  results,  then,  of  applying  the  new  crossflow  Reynolds 
number,  Eq.  (1),  to  these  profiles  are  also  found  in  Figure  10.  The  surface  of  Stetson's  cone  was 
cooled,  the  wall  of  King's  cone  was  adiabatic.  The  Stetson  data  were  taken  with  a  noisy 
free  stream. 

The  trends  of  the  three  sets  of  data  in  Figure  10  suggest  a  least  squares  fit  for  the 
correlation  of  R£.f<new)  and 

Quiet  (M=3.5) 

R,Knew)  =  26.7  38.0 

Noisy  (M=3.5) 

R,f(new)  =  14.7  +  33.9  W^^/Ue 
Noisy  (M=5.9) 

Ref(new)  =  26.5  +  25.2 

where  W^„AJe  is  in  %.  The  three  curves  are  plotted  in  Figure  1 1.  Below  the  curves,  the  flow 
has  not  undergone  transition.  The  values  of  R^Anew)  are  consistent  with  the  incompressible 
values  found  in  the  literature.  It  is  also  surprising  that  the  Stetson  data  agrees  so  well  with  the 
King  noisy  data. 

7.  CONCLUSIONS 

The  effects  of  nonequilibrium  chemistry  and  three-dimensionality  on  the  stability 
characteristics  of  hypersonic  and  supersonic  flows  have  been  discussed. 

In  2-D  and  axisymmetric  flows,  the  inclusion  of  chemistry  causes  a  shift  of  the  second 
mode  of  Mack  to  lower  frequencies.  This  is  found  to  be  due  to  the  increase  in  size  of  the  region 
of  relative  supersonic  flow  because  of  the  lower  speeds  of  sound  in  the  relatively  cooler 
boundary  layers. 

It  is  also  found  that  equilibrium  and  nonequilibrium  solutions  can  be  very  different 
depending  on  the  rates  of  the  reactions  relative  to  the  time  scales  of  convection  anu  diffusion.  In 
particular,  in  equilibrium-air  calculations,  modes  which  travel  supersonically  relative  to  the 
inviscid  region  are  shown  to  exist.  These  modes  are  a  superposition  of  incoming  and  outgoing 
disturbances  which  exhibit  oscillatory  behavior  because  of  the  finite  shock  stand-off  distance. 

For  the  examination  of  3-D  effects,  a  rotating  cone  is  successfully  used  as  a  model  of  a 
swept  wing.  A  large  parameter  range  of  flow  conditions  can  be  studied  and  trends  observed. 
The  value  of  this  should  become  more  apparent  from  the  subsequent  paragraphs. 

The  amplification  rate  of  the  first  mode  is  increased  by  a  factor  of  2  to  4  due  to  the 
presence  of  the  crossflow. 

The  most  unstable  crossflow  instability  has  a  nonzero  frequency. 

The  second  mode  in  a  3-D  boundary  layer  is  found  to  be  oblique  whereas  the  second 
mode  in  a  2-D  boundary  layer  is  2-D. 

Cooling  is  found  to  be  ineffective  in  controlling  crossflow. 

Considering  a  perfect  gas,  increasing  freestream  stagnation  temperature  only  slightly 
stabilizes  the  crossflow  instability. 

For  the  rotating  cone,  transition  location  (N  =  9)  does  not  correlate  with  the  traditional 
crossflow  Reynolds  number.  When  compressibility  and  cooling  effects  are  included,  Eq.  (1),  a 
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lelarionship  exists  between  the  new  crossflow  Reynolds  number  and  the  maximum  crossflow 
velocity  at  transition.  This  result  has  been  verified  with  the  yawed-cone  data  of  King  (1991)  and 
Stetson  (1982).  The  new  crossflow  Reynolds  number  is  calculated  solely  from  the  basic-state 
profiles  and,  as  such,  it  can  aid  in  transition  prediction  and  design  for  3-D  boundary  layers.  This 
formula  contains  no  geometry  explicitly  and  applies  for  flows  with  Prandtl  number  different 
from  unity. 
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9  FIGURES 


c)  Lociiions  of  gencr»]iz«i  infleciion  points  d)  TangcniitJ  velocity  at  generalized  inrieciion  points 


e)  Loeation  of  relative  sonic  points  for  neutral  disturbances 


f)  Neutral-disturbance  wavenumber  spacing  integrals 


Figure  1.  Basic-state  variations  with  respect  to  x 

•  •  -  Ideal  Gas  (Similarity  Solution),  - Ideal  Gas  (PNS), 

—  Equilibrium.  . Nonequilibrium 
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Figure  2.  Adiabatic  Hai-plaie  boundary.layer  ditturbajice-ttaie  ampUficaiion  rate*  ■  2000. 

Edge  Mach  number  e  10. 


Figure  A.  Routing  cone  in  supersonic  axial  flow  »nd  coordinate  system. 
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Fifure  7.  TradiiionaJ  crossnow  Reynolds  number  vs. 
shape  factor  ai  transition  for  various  wall-  and 
frcesircam-  temperature  condiliorts  at  M  =  3. 
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Figure  8.  Transition  Reynolds  number  vs.  cone 
rotational  speed  for  various  wall-  and  freesirerm- 
icmperature  condibons  at  M  =3. 
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Figure  9.  New  definition  of  crossflow  Reynolds  number 
(including  compressibility  effects)  versus  maximum 
crossHow  velocity  for  computational  rotating-cone 
data  assuming  a  transition  location  of  K=9.  Various 
frcesircam-  and  wall -temperature  conditions  as  well  as 
cone  geometries  arc  represented. 


(%) 

Figure  10.  New  definition  of  crossflow  Reynolds 
number  (  including  compressibility  effects  )  versus 
maximum  crossflow  velocity  for  experimental  yawed- 
cono  L.-ia  of  King  (1991)  at  M=3.5  and  Stetson  (1982) 
a;  M=i.9. 


(9i) 

Figure  11.  Linear  least-squares  approximation  of  new 
definition  of  crossflow  Reynolds  number  (including 
compressibility  effects)  as  a  function  of  maximum 
crossflow  vclocitv  for  experimental  vawed-cone  data  of 
King  (1991)  at  M=3.5  and  Stetson  (1982)  at  M=5.9. 
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10.  AVAILABLE  RESOURCES 

One  of  the  principal  strengths  of  our  team  at  Arizona  State  University  is  its  broad  skills  in 
analysis,  computations,  and  experiments.  We  facilitate  day-to-day  communication  between  the 
computational  work  and  the  experimental  work  through  two  IRIS  Graphics  Workstations  (3030 
and  3130)  and  two  DEC-5000  Workstations.  The  system,  with  state-of-the-art,  real-time,  three- 
dimensional,  color-graphics  software  (PLOT3D),  is  equipped  with  an  extensive  multi-user  and 
multi-task  environment  with  twelve  serial  lines.  Users  are  able  to  share  the  same  data  base  or 
experimental  information.  This  provides  the  hean  of  the  interaction  of  the  analytical, 
computational,  and  experimental  research. 

In  addition  to  the  super  computers  at  other  National  facilities,  the  network  includes 
access  to  the  IBM  4341A^M  and  Harris/VS  computers,  the  IBM  3090  Class  VI  machine,  and  the 
Cray  XMP  on  campus  as  well  as  the  MASSCOMP.  The  College  of  Engineering  at  ASU  is 
currently  also  equipped  with  several  VAX/780  and  VAX/785  minicomputers  exclusively  for 
research  purposes  (each  office  and  laboratory  has  a  hard-wired  RS232  interface).  These 
minicomputers  are  excellent  systems  for  program  development.  The  IRIS  and  DEC  machines 
can  access  all  the  features  available  in  those  minicomputers  through  the  existing  local  area 
networking  (Ethernet)  on  the  campus.  Furthermore,  the  system  can  communicate  directly  with 
other  research  facilities  to  share  information  through  telephone  couplings.  The  full  array  of 
computer  capabilities  from  super-mini  to  super-super  was  in  place  for  this  research. 
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